home *** CD-ROM | disk | FTP | other *** search
/ Enter 2004 January / enter-2004-01.iso / files / maxima-5.9.0.exe / {app} / share / maxima / 5.9.0 / src / numerical / slatec / djairy.lisp < prev    next >
Encoding:
Text File  |  2003-02-09  |  26.1 KB  |  495 lines

  1. ;;; Compiled by f2cl version 2.0 beta 2002-05-06
  2. ;;; 
  3. ;;; Options: ((:prune-labels nil) (:auto-save t) (:relaxed-array-decls t)
  4. ;;;           (:coerce-assigns :as-needed) (:array-type ':simple-array)
  5. ;;;           (:array-slicing nil) (:declare-common nil)
  6. ;;;           (:float-format double-float))
  7.  
  8. (in-package "SLATEC")
  9.  
  10.  
  11. (let ((n1 14)
  12.       (n2 23)
  13.       (n3 19)
  14.       (n4 15)
  15.       (m1 12)
  16.       (m2 21)
  17.       (m3 17)
  18.       (m4 13)
  19.       (fpi12 1.30899693899575)
  20.       (con2 5.03154716196777)
  21.       (con3 0.380004589867293)
  22.       (con4 0.833333333333333)
  23.       (con5 0.8660254037844392)
  24.       (ak1 (make-array 14 :element-type 'double-float))
  25.       (ak2 (make-array 23 :element-type 'double-float))
  26.       (ak3 (make-array 14 :element-type 'double-float))
  27.       (ajp (make-array 19 :element-type 'double-float))
  28.       (ajn (make-array 19 :element-type 'double-float))
  29.       (a (make-array 15 :element-type 'double-float))
  30.       (b (make-array 15 :element-type 'double-float))
  31.       (n1d 14)
  32.       (n2d 24)
  33.       (n3d 19)
  34.       (n4d 15)
  35.       (m1d 12)
  36.       (m2d 22)
  37.       (m3d 17)
  38.       (m4d 13)
  39.       (dak1 (make-array 14 :element-type 'double-float))
  40.       (dak2 (make-array 24 :element-type 'double-float))
  41.       (dak3 (make-array 14 :element-type 'double-float))
  42.       (dajp (make-array 19 :element-type 'double-float))
  43.       (dajn (make-array 19 :element-type 'double-float))
  44.       (da (make-array 15 :element-type 'double-float))
  45.       (db (make-array 15 :element-type 'double-float)))
  46.   (declare (type (simple-array double-float (24)) dak2)
  47.            (type (simple-array double-float (15)) db da b a)
  48.            (type (simple-array double-float (19)) dajn dajp ajn ajp)
  49.            (type (simple-array double-float (23)) ak2)
  50.            (type (simple-array double-float (14)) dak3 dak1 ak3 ak1)
  51.            (type double-float con5 con4 con3 con2 fpi12)
  52.            (type f2cl-lib:integer4 m4d m3d m2d m1d n4d n3d n2d n1d m4 m3 m2 m1
  53.             n4 n3 n2 n1))
  54.   (f2cl-lib:fset (f2cl-lib:fref ak1 (1) ((1 14))) 0.220423090987793)
  55.   (f2cl-lib:fset (f2cl-lib:fref ak1 (2) ((1 14))) -0.1252902427877)
  56.   (f2cl-lib:fset (f2cl-lib:fref ak1 (3) ((1 14))) 0.0103881163359194)
  57.   (f2cl-lib:fset (f2cl-lib:fref ak1 (4) ((1 14))) 8.228441520063431e-4)
  58.   (f2cl-lib:fset (f2cl-lib:fref ak1 (5) ((1 14))) -2.3461434589122598e-4)
  59.   (f2cl-lib:fset (f2cl-lib:fref ak1 (6) ((1 14))) 1.6382428017211603e-5)
  60.   (f2cl-lib:fset (f2cl-lib:fref ak1 (7) ((1 14))) 3.06902589573189e-7)
  61.   (f2cl-lib:fset (f2cl-lib:fref ak1 (8) ((1 14))) -1.29621999359332e-7)
  62.   (f2cl-lib:fset (f2cl-lib:fref ak1 (9) ((1 14))) 8.229081588236679e-9)
  63.   (f2cl-lib:fset (f2cl-lib:fref ak1 (10) ((1 14))) 1.5396396862329798e-11)
  64.   (f2cl-lib:fset (f2cl-lib:fref ak1 (11) ((1 14))) -3.3916546561568195e-11)
  65.   (f2cl-lib:fset (f2cl-lib:fref ak1 (12) ((1 14))) 2.0325325742362604e-12)
  66.   (f2cl-lib:fset (f2cl-lib:fref ak1 (13) ((1 14))) -1.10679546097884e-14)
  67.   (f2cl-lib:fset (f2cl-lib:fref ak1 (14) ((1 14))) -5.1616949778508e-15)
  68.   (f2cl-lib:fset (f2cl-lib:fref ak2 (1) ((1 23))) 0.274366150869598)
  69.   (f2cl-lib:fset (f2cl-lib:fref ak2 (2) ((1 23))) 0.005397909697369031)
  70.   (f2cl-lib:fset (f2cl-lib:fref ak2 (3) ((1 23))) -0.0015733922062118998)
  71.   (f2cl-lib:fset (f2cl-lib:fref ak2 (4) ((1 23))) 4.274275282487501e-4)
  72.   (f2cl-lib:fset (f2cl-lib:fref ak2 (5) ((1 23))) -1.12124917399925e-4)
  73.   (f2cl-lib:fset (f2cl-lib:fref ak2 (6) ((1 23))) 2.8876317131890405e-5)
  74.   (f2cl-lib:fset (f2cl-lib:fref ak2 (7) ((1 23))) -7.368042253705539e-6)
  75.   (f2cl-lib:fset (f2cl-lib:fref ak2 (8) ((1 23))) 1.87290209741024e-6)
  76.   (f2cl-lib:fset (f2cl-lib:fref ak2 (9) ((1 23))) -4.75892793962291e-7)
  77.   (f2cl-lib:fset (f2cl-lib:fref ak2 (10) ((1 23))) 1.2113041695590901e-7)
  78.   (f2cl-lib:fset (f2cl-lib:fref ak2 (11) ((1 23))) -3.09245374270614e-8)
  79.   (f2cl-lib:fset (f2cl-lib:fref ak2 (12) ((1 23))) 7.924547052826542e-9)
  80.   (f2cl-lib:fset (f2cl-lib:fref ak2 (13) ((1 23))) -2.0390244716791406e-9)
  81.   (f2cl-lib:fset (f2cl-lib:fref ak2 (14) ((1 23))) 5.268630565957421e-10)
  82.   (f2cl-lib:fset (f2cl-lib:fref ak2 (15) ((1 23))) -1.36704767639569e-10)
  83.   (f2cl-lib:fset (f2cl-lib:fref ak2 (16) ((1 23))) 3.561410390137079e-11)
  84.   (f2cl-lib:fset (f2cl-lib:fref ak2 (17) ((1 23))) -9.3138829654843e-12)
  85.   (f2cl-lib:fset (f2cl-lib:fref ak2 (18) ((1 23))) 2.44464450473635e-12)
  86.   (f2cl-lib:fset (f2cl-lib:fref ak2 (19) ((1 23))) -6.43840261990955e-13)
  87.   (f2cl-lib:fset (f2cl-lib:fref ak2 (20) ((1 23))) 1.7010603055934903e-13)
  88.   (f2cl-lib:fset (f2cl-lib:fref ak2 (21) ((1 23))) -4.5076010450328097e-14)
  89.   (f2cl-lib:fset (f2cl-lib:fref ak2 (22) ((1 23))) 1.19774799164811e-14)
  90.   (f2cl-lib:fset (f2cl-lib:fref ak2 (23) ((1 23))) -3.1907704086506605e-15)
  91.   (f2cl-lib:fset (f2cl-lib:fref ak3 (1) ((1 14))) 0.280271447340791)
  92.   (f2cl-lib:fset (f2cl-lib:fref ak3 (2) ((1 14))) -0.00178127042844379)
  93.   (f2cl-lib:fset (f2cl-lib:fref ak3 (3) ((1 14))) 4.034225796289991e-5)
  94.   (f2cl-lib:fset (f2cl-lib:fref ak3 (4) ((1 14))) -1.6324996526900298e-6)
  95.   (f2cl-lib:fset (f2cl-lib:fref ak3 (5) ((1 14))) 9.21181482476768e-8)
  96.   (f2cl-lib:fset (f2cl-lib:fref ak3 (6) ((1 14))) -6.522943302291551e-9)
  97.   (f2cl-lib:fset (f2cl-lib:fref ak3 (7) ((1 14))) 5.471384045765461e-10)
  98.   (f2cl-lib:fset (f2cl-lib:fref ak3 (8) ((1 14))) -5.244082518002599e-11)
  99.   (f2cl-lib:fset (f2cl-lib:fref ak3 (9) ((1 14))) 5.60477904117209e-12)
  100.   (f2cl-lib:fset (f2cl-lib:fref ak3 (10) ((1 14))) -6.56375244639313e-13)
  101.   (f2cl-lib:fset (f2cl-lib:fref ak3 (11) ((1 14))) 8.31285761966247e-14)
  102.   (f2cl-lib:fset (f2cl-lib:fref ak3 (12) ((1 14))) -1.12705134691063e-14)
  103.   (f2cl-lib:fset (f2cl-lib:fref ak3 (13) ((1 14))) 1.6226797659812903e-15)
  104.   (f2cl-lib:fset (f2cl-lib:fref ak3 (14) ((1 14))) -2.4648032431242603e-16)
  105.   (f2cl-lib:fset (f2cl-lib:fref ajp (1) ((1 19))) 0.0778952966437581)
  106.   (f2cl-lib:fset (f2cl-lib:fref ajp (2) ((1 19))) -0.184356363456801)
  107.   (f2cl-lib:fset (f2cl-lib:fref ajp (3) ((1 19))) 0.030141260521617398)
  108.   (f2cl-lib:fset (f2cl-lib:fref ajp (4) ((1 19))) 0.0305342724277608)
  109.   (f2cl-lib:fset (f2cl-lib:fref ajp (5) ((1 19))) -0.004954247025130791)
  110.   (f2cl-lib:fset (f2cl-lib:fref ajp (6) ((1 19))) -0.0017274955256395199)
  111.   (f2cl-lib:fset (f2cl-lib:fref ajp (7) ((1 19))) 2.4313763783918999e-4)
  112.   (f2cl-lib:fset (f2cl-lib:fref ajp (8) ((1 19))) 5.045647775170821e-5)
  113.   (f2cl-lib:fset (f2cl-lib:fref ajp (9) ((1 19))) -6.163165826952081e-6)
  114.   (f2cl-lib:fset (f2cl-lib:fref ajp (10) ((1 19))) -9.03986745510768e-7)
  115.   (f2cl-lib:fset (f2cl-lib:fref ajp (11) ((1 19))) 9.70243778355884e-8)
  116.   (f2cl-lib:fset (f2cl-lib:fref ajp (12) ((1 19))) 1.09639453305205e-8)
  117.   (f2cl-lib:fset (f2cl-lib:fref ajp (13) ((1 19))) -1.0471633058876599e-9)
  118.   (f2cl-lib:fset (f2cl-lib:fref ajp (14) ((1 19))) -9.603594413446459e-11)
  119.   (f2cl-lib:fset (f2cl-lib:fref ajp (15) ((1 19))) 8.253587894541341e-12)
  120.   (f2cl-lib:fset (f2cl-lib:fref ajp (16) ((1 19))) 6.36123439018768e-13)
  121.   (f2cl-lib:fset (f2cl-lib:fref ajp (17) ((1 19))) -4.9662961411601503e-14)
  122.   (f2cl-lib:fset (f2cl-lib:fref ajp (18) ((1 19))) -3.29810288929615e-15)
  123.   (f2cl-lib:fset (f2cl-lib:fref ajp (19) ((1 19))) 2.3579825203110397e-16)
  124.   (f2cl-lib:fset (f2cl-lib:fref ajn (1) ((1 19))) 0.038049788761724206)
  125.   (f2cl-lib:fset (f2cl-lib:fref ajn (2) ((1 19))) -0.24531954184554602)
  126.   (f2cl-lib:fset (f2cl-lib:fref ajn (3) ((1 19))) 0.165820623702696)
  127.   (f2cl-lib:fset (f2cl-lib:fref ajn (4) ((1 19))) 0.0749330045818789)
  128.   (f2cl-lib:fset (f2cl-lib:fref ajn (5) ((1 19))) -0.0263476288106641)
  129.   (f2cl-lib:fset (f2cl-lib:fref ajn (6) ((1 19))) -0.00592535597304981)
  130.   (f2cl-lib:fset (f2cl-lib:fref ajn (7) ((1 19))) 0.00144744409589804)
  131.   (f2cl-lib:fset (f2cl-lib:fref ajn (8) ((1 19))) 2.18311831322215e-4)
  132.   (f2cl-lib:fset (f2cl-lib:fref ajn (9) ((1 19))) -4.1066207768030405e-5)
  133.   (f2cl-lib:fset (f2cl-lib:fref ajn (10) ((1 19))) -4.6687499417176603e-6)
  134.   (f2cl-lib:fset (f2cl-lib:fref ajn (11) ((1 19))) 7.152188072771599e-7)
  135.   (f2cl-lib:fset (f2cl-lib:fref ajn (12) ((1 19))) 6.52964770854633e-8)
  136.   (f2cl-lib:fset (f2cl-lib:fref ajn (13) ((1 19))) -8.442840275659461e-9)
  137.   (f2cl-lib:fset (f2cl-lib:fref ajn (14) ((1 19))) -6.441861589769781e-10)
  138.   (f2cl-lib:fset (f2cl-lib:fref ajn (15) ((1 19))) 7.208022865052851e-11)
  139.   (f2cl-lib:fset (f2cl-lib:fref ajn (16) ((1 19))) 4.72465431717846e-12)
  140.   (f2cl-lib:fset (f2cl-lib:fref ajn (17) ((1 19))) -4.6602263254704496e-13)
  141.   (f2cl-lib:fset (f2cl-lib:fref ajn (18) ((1 19))) -2.6776271038918903e-14)
  142.   (f2cl-lib:fset (f2cl-lib:fref ajn (19) ((1 19))) 2.3616131657001907e-15)
  143.   (f2cl-lib:fset (f2cl-lib:fref a (1) ((1 15))) 0.490275424742791)
  144.   (f2cl-lib:fset (f2cl-lib:fref a (2) ((1 15))) 0.00157647277946204)
  145.   (f2cl-lib:fset (f2cl-lib:fref a (3) ((1 15))) -9.66195963140306e-5)
  146.   (f2cl-lib:fset (f2cl-lib:fref a (4) ((1 15))) 1.35916080268815e-7)
  147.   (f2cl-lib:fset (f2cl-lib:fref a (5) ((1 15))) 2.98157342654859e-7)
  148.   (f2cl-lib:fset (f2cl-lib:fref a (6) ((1 15))) -1.8682476755997898e-8)
  149.   (f2cl-lib:fset (f2cl-lib:fref a (7) ((1 15))) -1.0368573766714102e-9)
  150.   (f2cl-lib:fset (f2cl-lib:fref a (8) ((1 15))) 3.28660818434328e-10)
  151.   (f2cl-lib:fset (f2cl-lib:fref a (9) ((1 15))) -2.5709141063278e-11)
  152.   (f2cl-lib:fset (f2cl-lib:fref a (10) ((1 15))) -2.3235765530067695e-12)
  153.   (f2cl-lib:fset (f2cl-lib:fref a (11) ((1 15))) 9.575232790482551e-13)
  154.   (f2cl-lib:fset (f2cl-lib:fref a (12) ((1 15))) -1.2034082804971902e-13)
  155.   (f2cl-lib:fset (f2cl-lib:fref a (13) ((1 15))) -2.9090771677071503e-15)
  156.   (f2cl-lib:fset (f2cl-lib:fref a (14) ((1 15))) 4.5565645458014903e-15)
  157.   (f2cl-lib:fset (f2cl-lib:fref a (15) ((1 15))) -9.99003874810259e-16)
  158.   (f2cl-lib:fset (f2cl-lib:fref b (1) ((1 15))) 0.278593552803079)
  159.   (f2cl-lib:fset (f2cl-lib:fref b (2) ((1 15))) -0.00352915691882584)
  160.   (f2cl-lib:fset (f2cl-lib:fref b (3) ((1 15))) -2.31149677384994e-5)
  161.   (f2cl-lib:fset (f2cl-lib:fref b (4) ((1 15))) 4.7131784226356e-6)
  162.   (f2cl-lib:fset (f2cl-lib:fref b (5) ((1 15))) -1.12415907931333e-7)
  163.   (f2cl-lib:fset (f2cl-lib:fref b (6) ((1 15))) -2.0010030118433902e-8)
  164.   (f2cl-lib:fset (f2cl-lib:fref b (7) ((1 15))) 2.6094807530219305e-9)
  165.   (f2cl-lib:fset (f2cl-lib:fref b (8) ((1 15))) -3.55098136101216e-11)
  166.   (f2cl-lib:fset (f2cl-lib:fref b (9) ((1 15))) -3.50849978423875e-11)
  167.   (f2cl-lib:fset (f2cl-lib:fref b (10) ((1 15))) 5.830071879542019e-12)
  168.   (f2cl-lib:fset (f2cl-lib:fref b (11) ((1 15))) -2.04644828753326e-13)
  169.   (f2cl-lib:fset (f2cl-lib:fref b (12) ((1 15))) -1.1052917947674198e-13)
  170.   (f2cl-lib:fset (f2cl-lib:fref b (13) ((1 15))) 2.87724778038775e-14)
  171.   (f2cl-lib:fset (f2cl-lib:fref b (14) ((1 15))) -2.8820511100993906e-15)
  172.   (f2cl-lib:fset (f2cl-lib:fref b (15) ((1 15))) -3.32656311696166e-16)
  173.   (f2cl-lib:fset (f2cl-lib:fref dak1 (1) ((1 14))) 0.204567842307887)
  174.   (f2cl-lib:fset (f2cl-lib:fref dak1 (2) ((1 14))) -0.0661322739905664)
  175.   (f2cl-lib:fset (f2cl-lib:fref dak1 (3) ((1 14))) -0.00849845800989287)
  176.   (f2cl-lib:fset (f2cl-lib:fref dak1 (4) ((1 14))) 0.0031218349155628902)
  177.   (f2cl-lib:fset (f2cl-lib:fref dak1 (5) ((1 14))) -2.70016489829432e-4)
  178.   (f2cl-lib:fset (f2cl-lib:fref dak1 (6) ((1 14))) -6.356362986793869e-6)
  179.   (f2cl-lib:fset (f2cl-lib:fref dak1 (7) ((1 14))) 3.0239771240950897e-6)
  180.   (f2cl-lib:fset (f2cl-lib:fref dak1 (8) ((1 14))) -2.1831119533008797e-7)
  181.   (f2cl-lib:fset (f2cl-lib:fref dak1 (9) ((1 14))) -5.361942893328261e-10)
  182.   (f2cl-lib:fset (f2cl-lib:fref dak1 (10) ((1 14))) 1.1309803562231e-9)
  183.   (f2cl-lib:fset (f2cl-lib:fref dak1 (11) ((1 14))) -7.43023834629073e-11)
  184.   (f2cl-lib:fset (f2cl-lib:fref dak1 (12) ((1 14))) 4.2880417082689104e-13)
  185.   (f2cl-lib:fset (f2cl-lib:fref dak1 (13) ((1 14))) 2.23810925754539e-13)
  186.   (f2cl-lib:fset (f2cl-lib:fref dak1 (14) ((1 14))) -1.39140135641182e-14)
  187.   (f2cl-lib:fset (f2cl-lib:fref dak2 (1) ((1 24))) 0.29333234388323)
  188.   (f2cl-lib:fset (f2cl-lib:fref dak2 (2) ((1 24))) -0.00806196784743112)
  189.   (f2cl-lib:fset (f2cl-lib:fref dak2 (3) ((1 24))) 0.0024254017233314)
  190.   (f2cl-lib:fset (f2cl-lib:fref dak2 (4) ((1 24))) -6.82297548850235e-4)
  191.   (f2cl-lib:fset (f2cl-lib:fref dak2 (5) ((1 24))) 1.85786427751181e-4)
  192.   (f2cl-lib:fset (f2cl-lib:fref dak2 (6) ((1 24))) -4.974574476840591e-5)
  193.   (f2cl-lib:fset (f2cl-lib:fref dak2 (7) ((1 24))) 1.3209068123949702e-5)
  194.   (f2cl-lib:fset (f2cl-lib:fref dak2 (8) ((1 24))) -3.4952824044494296e-6)
  195.   (f2cl-lib:fset (f2cl-lib:fref dak2 (9) ((1 24))) 9.24362451078835e-7)
  196.   (f2cl-lib:fset (f2cl-lib:fref dak2 (10) ((1 24))) -2.4473267152186698e-7)
  197.   (f2cl-lib:fset (f2cl-lib:fref dak2 (11) ((1 24))) 6.4930783764891e-8)
  198.   (f2cl-lib:fset (f2cl-lib:fref dak2 (12) ((1 24))) -1.72717621501538e-8)
  199.   (f2cl-lib:fset (f2cl-lib:fref dak2 (13) ((1 24))) 4.6072576360465606e-9)
  200.   (f2cl-lib:fset (f2cl-lib:fref dak2 (14) ((1 24))) -1.2324905529154999e-9)
  201.   (f2cl-lib:fset (f2cl-lib:fref dak2 (15) ((1 24))) 3.30620409488102e-10)
  202.   (f2cl-lib:fset (f2cl-lib:fref dak2 (16) ((1 24))) -8.892520997724009e-11)
  203.   (f2cl-lib:fset (f2cl-lib:fref dak2 (17) ((1 24))) 2.3977331987829795e-11)
  204.   (f2cl-lib:fset (f2cl-lib:fref dak2 (18) ((1 24))) -6.480139211534499e-12)
  205.   (f2cl-lib:fset (f2cl-lib:fref dak2 (19) ((1 24))) 1.75510132023731e-12)
  206.   (f2cl-lib:fset (f2cl-lib:fref dak2 (20) ((1 24))) -4.7630382983363695e-13)
  207.   (f2cl-lib:fset (f2cl-lib:fref dak2 (21) ((1 24))) 1.2949824110081e-13)
  208.   (f2cl-lib:fset (f2cl-lib:fref dak2 (22) ((1 24))) -3.5267962221043e-14)
  209.   (f2cl-lib:fset (f2cl-lib:fref dak2 (23) ((1 24))) 9.62005151585923e-15)
  210.   (f2cl-lib:fset (f2cl-lib:fref dak2 (24) ((1 24))) -2.62786914342292e-15)
  211.   (f2cl-lib:fset (f2cl-lib:fref dak3 (1) ((1 14))) 0.284675828811349)
  212.   (f2cl-lib:fset (f2cl-lib:fref dak3 (2) ((1 14))) 0.0025307307261908003)
  213.   (f2cl-lib:fset (f2cl-lib:fref dak3 (3) ((1 14))) -4.83481130337976e-5)
  214.   (f2cl-lib:fset (f2cl-lib:fref dak3 (4) ((1 14))) 1.8490728394634298e-6)
  215.   (f2cl-lib:fset (f2cl-lib:fref dak3 (5) ((1 14))) -1.01418491178576e-7)
  216.   (f2cl-lib:fset (f2cl-lib:fref dak3 (6) ((1 14))) 7.059256344571531e-9)
  217.   (f2cl-lib:fset (f2cl-lib:fref dak3 (7) ((1 14))) -5.85325291400382e-10)
  218.   (f2cl-lib:fset (f2cl-lib:fref dak3 (8) ((1 14))) 5.563576888313391e-11)
  219.   (f2cl-lib:fset (f2cl-lib:fref dak3 (9) ((1 14))) -5.908890947794999e-12)
  220.   (f2cl-lib:fset (f2cl-lib:fref dak3 (10) ((1 14))) 6.88574353784436e-13)
  221.   (f2cl-lib:fset (f2cl-lib:fref dak3 (11) ((1 14))) -8.68588256452194e-14)
  222.   (f2cl-lib:fset (f2cl-lib:fref dak3 (12) ((1 14))) 1.17374762617213e-14)
  223.   (f2cl-lib:fset (f2cl-lib:fref dak3 (13) ((1 14))) -1.68523146510923e-15)
  224.   (f2cl-lib:fset (f2cl-lib:fref dak3 (14) ((1 14))) 2.5537477309705603e-16)
  225.   (f2cl-lib:fset (f2cl-lib:fref dajp (1) ((1 19))) 0.0653219131311457)
  226.   (f2cl-lib:fset (f2cl-lib:fref dajp (2) ((1 19))) -0.120262933688823)
  227.   (f2cl-lib:fset (f2cl-lib:fref dajp (3) ((1 19))) 0.00978010236263823)
  228.   (f2cl-lib:fset (f2cl-lib:fref dajp (4) ((1 19))) 0.016794842923050498)
  229.   (f2cl-lib:fset (f2cl-lib:fref dajp (5) ((1 19))) -0.00197146140182132)
  230.   (f2cl-lib:fset (f2cl-lib:fref dajp (6) ((1 19))) -8.45560295098867e-4)
  231.   (f2cl-lib:fset (f2cl-lib:fref dajp (7) ((1 19))) 9.428896207019759e-5)
  232.   (f2cl-lib:fset (f2cl-lib:fref dajp (8) ((1 19))) 2.25827860945475e-5)
  233.   (f2cl-lib:fset (f2cl-lib:fref dajp (9) ((1 19))) -2.29067870915987e-6)
  234.   (f2cl-lib:fset (f2cl-lib:fref dajp (10) ((1 19))) -3.7634399113691896e-7)
  235.   (f2cl-lib:fset (f2cl-lib:fref dajp (11) ((1 19))) 3.4566393355956504e-8)
  236.   (f2cl-lib:fset (f2cl-lib:fref dajp (12) ((1 19))) 4.2961133200300705e-9)
  237.   (f2cl-lib:fset (f2cl-lib:fref dajp (13) ((1 19))) -3.58673691214989e-10)
  238.   (f2cl-lib:fset (f2cl-lib:fref dajp (14) ((1 19))) -3.5724588136189495e-11)
  239.   (f2cl-lib:fset (f2cl-lib:fref dajp (15) ((1 19))) 2.7269609106633597e-12)
  240.   (f2cl-lib:fset (f2cl-lib:fref dajp (16) ((1 19))) 2.2612065309577098e-13)
  241.   (f2cl-lib:fset (f2cl-lib:fref dajp (17) ((1 19))) -1.5876320523830303e-14)
  242.   (f2cl-lib:fset (f2cl-lib:fref dajp (18) ((1 19))) -1.1260437448512499e-15)
  243.   (f2cl-lib:fset (f2cl-lib:fref dajp (19) ((1 19))) 7.313275295153671e-17)
  244.   (f2cl-lib:fset (f2cl-lib:fref dajn (1) ((1 19))) 0.0108594539632967)
  245.   (f2cl-lib:fset (f2cl-lib:fref dajn (2) ((1 19))) 0.0853313194857091)
  246.   (f2cl-lib:fset (f2cl-lib:fref dajn (3) ((1 19))) -0.315277068113058)
  247.   (f2cl-lib:fset (f2cl-lib:fref dajn (4) ((1 19))) -0.0878420725294257)
  248.   (f2cl-lib:fset (f2cl-lib:fref dajn (5) ((1 19))) 0.0553251906976048)
  249.   (f2cl-lib:fset (f2cl-lib:fref dajn (6) ((1 19))) 0.009416740605032411)
  250.   (f2cl-lib:fset (f2cl-lib:fref dajn (7) ((1 19))) -0.00332187026018996)
  251.   (f2cl-lib:fset (f2cl-lib:fref dajn (8) ((1 19))) -4.1115734315682606e-4)
  252.   (f2cl-lib:fset (f2cl-lib:fref dajn (9) ((1 19))) 1.0129732689134603e-4)
  253.   (f2cl-lib:fset (f2cl-lib:fref dajn (10) ((1 19))) 9.876336822083961e-6)
  254.   (f2cl-lib:fset (f2cl-lib:fref dajn (11) ((1 19))) -1.8731296981239298e-6)
  255.   (f2cl-lib:fset (f2cl-lib:fref dajn (12) ((1 19))) -1.5079850013146798e-7)
  256.   (f2cl-lib:fset (f2cl-lib:fref dajn (13) ((1 19))) 2.3268766952539402e-8)
  257.   (f2cl-lib:fset (f2cl-lib:fref dajn (14) ((1 19))) 1.5959991741922502e-9)
  258.   (f2cl-lib:fset (f2cl-lib:fref dajn (15) ((1 19))) -2.0766592266838502e-10)
  259.   (f2cl-lib:fset (f2cl-lib:fref dajn (16) ((1 19))) -1.24103350500302e-11)
  260.   (f2cl-lib:fset (f2cl-lib:fref dajn (17) ((1 19))) 1.39631765331043e-12)
  261.   (f2cl-lib:fset (f2cl-lib:fref dajn (18) ((1 19))) 7.3940097115574e-14)
  262.   (f2cl-lib:fset (f2cl-lib:fref dajn (19) ((1 19))) -7.328874756275e-15)
  263.   (f2cl-lib:fset (f2cl-lib:fref da (1) ((1 15))) 0.49162732110460106)
  264.   (f2cl-lib:fset (f2cl-lib:fref da (2) ((1 15))) 0.00311164930427489)
  265.   (f2cl-lib:fset (f2cl-lib:fref da (3) ((1 15))) 8.23140762854081e-5)
  266.   (f2cl-lib:fset (f2cl-lib:fref da (4) ((1 15))) -4.6176977617214193e-6)
  267.   (f2cl-lib:fset (f2cl-lib:fref da (5) ((1 15))) -6.13158880534626e-8)
  268.   (f2cl-lib:fset (f2cl-lib:fref da (6) ((1 15))) 2.8729580465652e-8)
  269.   (f2cl-lib:fset (f2cl-lib:fref da (7) ((1 15))) -1.8195971537211703e-9)
  270.   (f2cl-lib:fset (f2cl-lib:fref da (8) ((1 15))) -1.44752826642035e-10)
  271.   (f2cl-lib:fset (f2cl-lib:fref da (9) ((1 15))) 4.5372404342042194e-11)
  272.   (f2cl-lib:fset (f2cl-lib:fref da (10) ((1 15))) -3.99655065847223e-12)
  273.   (f2cl-lib:fset (f2cl-lib:fref da (11) ((1 15))) -3.24089119830323e-13)
  274.   (f2cl-lib:fset (f2cl-lib:fref da (12) ((1 15))) 1.6209895256874102e-13)
  275.   (f2cl-lib:fset (f2cl-lib:fref da (13) ((1 15))) -2.4076524797405702e-14)
  276.   (f2cl-lib:fset (f2cl-lib:fref da (14) ((1 15))) 1.6938481128449098e-16)
  277.   (f2cl-lib:fset (f2cl-lib:fref da (15) ((1 15))) 8.17900786477396e-16)
  278.   (f2cl-lib:fset (f2cl-lib:fref db (1) ((1 15))) -0.277571356944231)
  279.   (f2cl-lib:fset (f2cl-lib:fref db (2) ((1 15))) 0.0044421283341992005)
  280.   (f2cl-lib:fset (f2cl-lib:fref db (3) ((1 15))) -8.42328522190089e-5)
  281.   (f2cl-lib:fset (f2cl-lib:fref db (4) ((1 15))) -2.5804031841870995e-6)
  282.   (f2cl-lib:fset (f2cl-lib:fref db (5) ((1 15))) 3.4238972021762104e-7)
  283.   (f2cl-lib:fset (f2cl-lib:fref db (6) ((1 15))) -6.24286894709776e-9)
  284.   (f2cl-lib:fset (f2cl-lib:fref db (7) ((1 15))) -2.3637783684457703e-9)
  285.   (f2cl-lib:fset (f2cl-lib:fref db (8) ((1 15))) 3.16991042656673e-10)
  286.   (f2cl-lib:fset (f2cl-lib:fref db (9) ((1 15))) -4.4099569165819097e-12)
  287.   (f2cl-lib:fset (f2cl-lib:fref db (10) ((1 15))) -5.186742210935751e-12)
  288.   (f2cl-lib:fset (f2cl-lib:fref db (11) ((1 15))) 9.64874015137022e-13)
  289.   (f2cl-lib:fset (f2cl-lib:fref db (12) ((1 15))) -4.9019057660870996e-14)
  290.   (f2cl-lib:fset (f2cl-lib:fref db (13) ((1 15))) -1.77253430678112e-14)
  291.   (f2cl-lib:fset (f2cl-lib:fref db (14) ((1 15))) 5.55950610442662e-15)
  292.   (f2cl-lib:fset (f2cl-lib:fref db (15) ((1 15))) -7.117933375795299e-16)
  293.   (defun djairy (x rx c ai dai)
  294.     (declare (type double-float dai ai c rx x))
  295.     (prog ((ccv 0.0) (cv 0.0) (ec 0.0) (e1 0.0) (e2 0.0) (f1 0.0) (f2 0.0)
  296.            (rtrx 0.0) (scv 0.0) (t_ 0.0) (temp1 0.0) (temp2 0.0) (tt 0.0) (i 0)
  297.            (j 0))
  298.       (declare (type f2cl-lib:integer4 j i)
  299.                (type double-float tt temp2 temp1 t_ scv rtrx f2 f1 e2 e1 ec cv
  300.                 ccv))
  301.       (if (< x 0.0) (go label90))
  302.       (if (> c 5.0) (go label60))
  303.       (if (> x 1.2) (go label30))
  304.       (setf t_ (* (- (+ x x) 1.2) con4))
  305.       (setf tt (+ t_ t_))
  306.       (setf j n1)
  307.       (setf f1 (f2cl-lib:fref ak1 (j) ((1 14))))
  308.       (setf f2 0.0)
  309.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  310.                     ((> i m1) nil)
  311.         (tagbody
  312.           (setf j (f2cl-lib:int-sub j 1))
  313.           (setf temp1 f1)
  314.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref ak1 (j) ((1 14)))))
  315.           (setf f2 temp1)
  316.          label10))
  317.       (setf ai (+ (- (* t_ f1) f2) (f2cl-lib:fref ak1 (1) ((1 14)))))
  318.       (setf j n1d)
  319.       (setf f1 (f2cl-lib:fref dak1 (j) ((1 14))))
  320.       (setf f2 0.0)
  321.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  322.                     ((> i m1d) nil)
  323.         (tagbody
  324.           (setf j (f2cl-lib:int-sub j 1))
  325.           (setf temp1 f1)
  326.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dak1 (j) ((1 14)))))
  327.           (setf f2 temp1)
  328.          label20))
  329.       (setf dai (- (+ (- (* t_ f1) f2) (f2cl-lib:fref dak1 (1) ((1 14))))))
  330.       (go end_label)
  331.      label30
  332.       (setf t_ (* (- (+ x x) con2) con3))
  333.       (setf tt (+ t_ t_))
  334.       (setf j n2)
  335.       (setf f1 (f2cl-lib:fref ak2 (j) ((1 23))))
  336.       (setf f2 0.0)
  337.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  338.                     ((> i m2) nil)
  339.         (tagbody
  340.           (setf j (f2cl-lib:int-sub j 1))
  341.           (setf temp1 f1)
  342.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref ak2 (j) ((1 23)))))
  343.           (setf f2 temp1)
  344.          label40))
  345.       (setf rtrx (f2cl-lib:fsqrt rx))
  346.       (setf ec (exp (- c)))
  347.       (setf ai
  348.               (/ (* ec (+ (- (* t_ f1) f2) (f2cl-lib:fref ak2 (1) ((1 23)))))
  349.                  rtrx))
  350.       (setf j n2d)
  351.       (setf f1 (f2cl-lib:fref dak2 (j) ((1 24))))
  352.       (setf f2 0.0)
  353.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  354.                     ((> i m2d) nil)
  355.         (tagbody
  356.           (setf j (f2cl-lib:int-sub j 1))
  357.           (setf temp1 f1)
  358.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dak2 (j) ((1 24)))))
  359.           (setf f2 temp1)
  360.          label50))
  361.       (setf dai
  362.               (* (- ec)
  363.                  (+ (- (* t_ f1) f2) (f2cl-lib:fref dak2 (1) ((1 24))))
  364.                  rtrx))
  365.       (go end_label)
  366.      label60
  367.       (setf t_ (- (/ 10.0 c) 1.0))
  368.       (setf tt (+ t_ t_))
  369.       (setf j n1)
  370.       (setf f1 (f2cl-lib:fref ak3 (j) ((1 14))))
  371.       (setf f2 0.0)
  372.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  373.                     ((> i m1) nil)
  374.         (tagbody
  375.           (setf j (f2cl-lib:int-sub j 1))
  376.           (setf temp1 f1)
  377.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref ak3 (j) ((1 14)))))
  378.           (setf f2 temp1)
  379.          label70))
  380.       (setf rtrx (f2cl-lib:fsqrt rx))
  381.       (setf ec (exp (- c)))
  382.       (setf ai
  383.               (/ (* ec (+ (- (* t_ f1) f2) (f2cl-lib:fref ak3 (1) ((1 14)))))
  384.                  rtrx))
  385.       (setf j n1d)
  386.       (setf f1 (f2cl-lib:fref dak3 (j) ((1 14))))
  387.       (setf f2 0.0)
  388.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  389.                     ((> i m1d) nil)
  390.         (tagbody
  391.           (setf j (f2cl-lib:int-sub j 1))
  392.           (setf temp1 f1)
  393.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dak3 (j) ((1 14)))))
  394.           (setf f2 temp1)
  395.          label80))
  396.       (setf dai
  397.               (* (- rtrx)
  398.                  ec
  399.                  (+ (- (* t_ f1) f2) (f2cl-lib:fref dak3 (1) ((1 14))))))
  400.       (go end_label)
  401.      label90
  402.       (if (> c 5.0) (go label120))
  403.       (setf t_ (- (* 0.4 c) 1.0))
  404.       (setf tt (+ t_ t_))
  405.       (setf j n3)
  406.       (setf f1 (f2cl-lib:fref ajp (j) ((1 19))))
  407.       (setf e1 (f2cl-lib:fref ajn (j) ((1 19))))
  408.       (setf f2 0.0)
  409.       (setf e2 0.0)
  410.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  411.                     ((> i m3) nil)
  412.         (tagbody
  413.           (setf j (f2cl-lib:int-sub j 1))
  414.           (setf temp1 f1)
  415.           (setf temp2 e1)
  416.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref ajp (j) ((1 19)))))
  417.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref ajn (j) ((1 19)))))
  418.           (setf f2 temp1)
  419.           (setf e2 temp2)
  420.          label100))
  421.       (setf ai
  422.               (- (+ (- (* t_ e1) e2) (f2cl-lib:fref ajn (1) ((1 19))))
  423.                  (* x (+ (- (* t_ f1) f2) (f2cl-lib:fref ajp (1) ((1 19)))))))
  424.       (setf j n3d)
  425.       (setf f1 (f2cl-lib:fref dajp (j) ((1 19))))
  426.       (setf e1 (f2cl-lib:fref dajn (j) ((1 19))))
  427.       (setf f2 0.0)
  428.       (setf e2 0.0)
  429.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  430.                     ((> i m3d) nil)
  431.         (tagbody
  432.           (setf j (f2cl-lib:int-sub j 1))
  433.           (setf temp1 f1)
  434.           (setf temp2 e1)
  435.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref dajp (j) ((1 19)))))
  436.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref dajn (j) ((1 19)))))
  437.           (setf f2 temp1)
  438.           (setf e2 temp2)
  439.          label110))
  440.       (setf dai
  441.               (+ (* x x (+ (- (* t_ f1) f2) (f2cl-lib:fref dajp (1) ((1 19)))))
  442.                  (+ (- (* t_ e1) e2) (f2cl-lib:fref dajn (1) ((1 19))))))
  443.       (go end_label)
  444.      label120
  445.       (setf t_ (- (/ 10.0 c) 1.0))
  446.       (setf tt (+ t_ t_))
  447.       (setf j n4)
  448.       (setf f1 (f2cl-lib:fref a (j) ((1 15))))
  449.       (setf e1 (f2cl-lib:fref b (j) ((1 15))))
  450.       (setf f2 0.0)
  451.       (setf e2 0.0)
  452.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  453.                     ((> i m4) nil)
  454.         (tagbody
  455.           (setf j (f2cl-lib:int-sub j 1))
  456.           (setf temp1 f1)
  457.           (setf temp2 e1)
  458.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref a (j) ((1 15)))))
  459.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref b (j) ((1 15)))))
  460.           (setf f2 temp1)
  461.           (setf e2 temp2)
  462.          label130))
  463.       (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref a (1) ((1 15)))))
  464.       (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref b (1) ((1 15)))))
  465.       (setf rtrx (f2cl-lib:fsqrt rx))
  466.       (setf cv (- c fpi12))
  467.       (setf ccv (cos cv))
  468.       (setf scv (sin cv))
  469.       (setf ai (/ (- (* temp1 ccv) (* temp2 scv)) rtrx))
  470.       (setf j n4d)
  471.       (setf f1 (f2cl-lib:fref da (j) ((1 15))))
  472.       (setf e1 (f2cl-lib:fref db (j) ((1 15))))
  473.       (setf f2 0.0)
  474.       (setf e2 0.0)
  475.       (f2cl-lib:fdo (i 1 (f2cl-lib:int-add i 1))
  476.                     ((> i m4d) nil)
  477.         (tagbody
  478.           (setf j (f2cl-lib:int-sub j 1))
  479.           (setf temp1 f1)
  480.           (setf temp2 e1)
  481.           (setf f1 (+ (- (* tt f1) f2) (f2cl-lib:fref da (j) ((1 15)))))
  482.           (setf e1 (+ (- (* tt e1) e2) (f2cl-lib:fref db (j) ((1 15)))))
  483.           (setf f2 temp1)
  484.           (setf e2 temp2)
  485.          label140))
  486.       (setf temp1 (+ (- (* t_ f1) f2) (f2cl-lib:fref da (1) ((1 15)))))
  487.       (setf temp2 (+ (- (* t_ e1) e2) (f2cl-lib:fref db (1) ((1 15)))))
  488.       (setf e1 (+ (* ccv con5) (* 0.5 scv)))
  489.       (setf e2 (- (* scv con5) (* 0.5 ccv)))
  490.       (setf dai (* (- (* temp1 e1) (* temp2 e2)) rtrx))
  491.       (go end_label)
  492.      end_label
  493.       (return (values nil nil nil ai dai)))))
  494.  
  495.